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Abstract 

Continuing work initiated in an earlier publication [Torigoe et al. Phys. Rev. Lett. 102, 
251101 (2009)], gravitational wave forms for a three-body system in Lagrange's orbit are considered 
especially in an analytic method. First, we derive an expression of the three-body wave forms 
at the mass quadrupole, octupole and current quadrupole orders. By using the expressions, we 
solve a gravitational-wave inverse problem of determining the source parameters to this particular 
configuration (three masses, a distance of the source to an observer, and the orbital inclination 
angle to the line of sight) through observations of the gravitational wave forms alone. For this 
purpose, the chirp mass to a three-body system in the particular configuration is expressed in 
terms of only the mass ratios by deleting initial angle positions. We discuss also whether and how 
a binary source can be distinguished from a three-body system in Lagrange's orbit or others. 
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I. INTRODUCTION 



"Can one hear the shape of a drum? 1 is a famous question posed by Mark Kac in 1966 
where to "hear" the shape of a drum is to infer information about the shape of the 
drumhead from the sound it makes. Such a question can be traced back to Hermann Weyl 



Now, it is interesting to pose a gravitational-wave inverse problem 



? or the forthcoming 



il9|. To 
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gravitational-wave astronomy by ground-based or space-borne detectors 
"hear" a source through gravitational-wave observations is to extract the information about 
the source from the gravitational waves it makes. 

It is of general interest to ask "can one tell how many apples are falling in the dark of 
night? 1 One simpler question is how and whether two-body and three-body gravitating 
systems can be distinguished through observations of gravitational waves that are made 
by these sources. Recently, this issue has been addressed [h]]. They found that there is 
a case that quadrupole wave forms for two-body and three-body gravitating systems (in 
a particular configuration) cannot be distinguished even with observing the chirp, though 
different numbers of self-gravitating particles (in different types of periodic motion) are 
most likely to generate very different shapes of gravitational waves. In order to break this 
degeneracy between two-body and three-body sources in the particular configuration for 
the quadrupolar wave (even with observing frequency sweep), they suggested that higher 
multipolar wave forms (octupole for their cases) are needed. 

The purpose of the present paper is to investigate gravitational wave forms for a three- 
body system in the particular configuration (as Lagrange's orbit) in more detail, especially in 
analytic manners up to the mass octupole and current quadrupole orders. In particular, we 
shall discuss whether and how the source parameters (three masses, a distance of the sources 
to an observer, and their orbital inclination angle to the line of sight) can be determined 
through observations of the gravitational waveforms alone. 

Inspiraling and finally merging binary compact stars are most likely astrophysical sources 
for the forthcoming direct detections of gravitational ripples (and consequently gravitational 

n n n u n 

waves astronomy) by a lot of efforts by the on-going or designed detectors [4., |5|, |6|, |7J, la, 
M. Merging neutron stars and black holes have been successfully simulated by numerical 

^nnn 

131 . 114 , 1 1 51 ] - Furthermore, analytic methods also have provided accurate 



relativity [11 
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waveform templates for inspiraling compact binaries, by post-Newtonian approaches (See 
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IS Il7] for reviews) and also by black hole perturbation techniques especially at the linear 
order in mass ratio (See also 18] for reviews). Bridges between the inspiraling stage and 
final merging phase are currently under construction (e.g., HQ). 

There is a growing interest in pot ential astrophysical sources of gravitational waves in- 
volving 3-body interactions (e.g., 1(3, 21, 22, 23] and references therein). Even the classical 



three-body (or N-body) problem in Newtonian gravity admits an infinite number of so- 
lutions, some of which express regular orbits and the others are chaotic, and in fact an 
increasing number of periodic orbits are found 2J, |25j. For the sake of simplicity, we focus 



on Lagrange's equilateral triangle solution, mostly because it offers a nice model tractable 
completely by hand (See below for more details). 



The linear perturbation analysis in Newton gravity [24| shows that Lagrange triangular 
points L4 and L5 to the restricted three-body problem, in which one of three bodies is 
assumed as a test mass (e.g., an asteroid in the solar system), is stable if the mass ratio of 
the remaining two bodies is less than 0.0385, though relativistic corrections to stability of 
the system are poorly understood. Indeed, the ratio of the Jovian mass to the solar mass is 
O(10~ 3 ). 

Here, it is worthwhile to mention previous works. Nakamura and Oohara [26] studied 
numerically the luminosity of gravitational radiation by N test particles orbiting around a 
Schwarzschild black hole, as an extension of Detweiler's analysis of the N — 1 case [^tJ by 
using Teukolsky equation 28j, in order to show the phase cancellation effect, which had been 
pointed out by Nakamura and Sasaki [29( . It should be noted that their N particles are test 
masses but not self-gravitating. The aim and setting of the present paper are different from 
those previous works. 

In this paper, we shall present an expression of gravitational waves for a three-body 
system in Lagrange's orbit, especially as a function of mass ratios. In addition, we shall 
discuss how to determine the source parameters (three masses, a distance of the source to 
an observer, and the orbital inclination angle to the line of sight) through observations of 
the gravitational waveforms alone. We also discuss a possible test; which source is a binary, 
a three-body system in Lagrange's orbit or others? 

This paper is organized as follows. In section 2, we shall briefly summarize a notation 
and formulation for describing a three-body system in Lagrange's orbit and its gravitational 
waves. Wave forms at the mass quadrupole, mass octupole and current quadrupole orders 
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are obtained in an analytic method. In section 3, we discuss a method of determining the 
source parameters from gravitational- wave observations alone. A binary source test is also 
discussed. Section 4 is devoted to the conclusion. In Appendix, we shall present calculations 
for the chirp mass and some useful relations. Throughout this paper, we take the units of 
G = c = 1. Latin indices take 1, 2, 3, except for p which labels each body. 



II. NOTATION AND FORMULATION 



A. Linearized gravitational waves 



For a wave propagation direction denoted by a unit vector n a , we define the transverse- 
traceless projection operator as 



P b = 5" - n rf 

1 a u a "'a.' 1 i 



(1) 



where 5 b a denotes the Kronecker's delta symbol. The linearized waves in the wave zone 
are expressed in terms of mass and current multipole moments denoted by I At and S Ae , 

ri 

respectively [30j. These radiative multipole moments at the Newtonian order are related 
with the source position x a , mass density p and velocity v a as 

- STF 

pX Ae d 3 x 



S A * 



e a l c x b pv c X A ^d 3 x 



STF 



(2) 
(3) 



where STF denotes the symmetric tracefree part, e a b c denotes the Levi-Civita symbol in 
a three-dimensional Euclidean space, and we define the product of i spatial coordinates as 
X Ai = x ai x a2 ■ ■ -x ai . We consider a system of spherical bodies approximated by massive 
particles. Then, Eqs. (J2J) and (jHJ) become 

STF 

(4) 



S A * 



N 



P =i 

N 



P =i 



STF 



(5) 



where N denotes the number of the particles and the subscript p denotes the p-th body. 
In terms of the multipole moments, the linearized waves at the wave zone are expressed 
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as 3 



hjT(t,x) 



£ 

OO 



,5 



(t - r)iV^- 2 



TT 



(6) 



where t and r mean time and source distance, respectively, in the Minkowskian spher- 
ical coordinates (t, r, 6,$), TT denotes the transverse-traceless part, (£) denotes the t- 
th time derivative, and we define the tensor product of t — 2 unit radial vectors as 

= n ai n a.2 . . . n a, e _ 2 _ 

For a binary case, mass quadrupolar, octupolar and current quadrupolar waves were 
considered fully by Blanchet and Schafer 



3l|, where they showed that the octupolar part is 



linearly proportional to a mass difference. 



B. Lagrange's solution 

Let us consider the Lagrange's solution for a three-body system (on x-y plane), where 
each mass is denoted by m p (p = 1,2,3). The initial positions of each mass are expressed 
as Xi = (0, 0) , = a(v3/2, 1/2), and x% = a(0, 1), where the side of a regular triangle is 
denoted as a [lol. . 

We choose the spatial coordinates such that the center of mass (COM) is at rest as 
(^com,2/com) = a(y/3u 2 /2, {y 2 + v^)/2)), where the total mass and mass ratio are denoted 
as m tot = J2 P m p an d v p = m P / m tot, respectively. We have an identity as 

v\ + v 2 + v 3 = 1. (7) 

The orbital frequency u for the triangle satisfies 

w = —3-, (8) 

which takes the same form as Kepler's third law for a binary system but with the total mass 
of three masses. 

Henceforth, it is convenient to employ the COM coordinates (X, Y) that can be obtained 
by a translation from (x, y). In the COM coordinates, the location of each mass is expressed 
as 

X p = a p (cos(ut + 9 p ),sm(ujt + 9 P )), (9) 



m 2 ». y 
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FIG. 1: Definition of 9 P in the Lagrange's orbit. The angle 9 P is measured from X-axis to the 
direction of each mass at the initial time. 

where a p is defined as oi = V^com + Vcom, °2 = V (3 1/2 a/2 - i C om) 2 + (a/2 - ?/com) 2 , 
and 03 = a/ x com + (a — 2/com) 2 , respectively, and 9 P denotes the angle between the new X- 
axis and the direction of each mass at t — (See Fig. [T]). In practice, computations can be 
simplified by using complex variables. In particular, variables for the triangle configuration 
8 p are written in terms of the mass ratios (See Appendix for more detail). 



III. WAVE FORMS FOR THREE BODIES IN LAGRANGE'S ORBIT 



Let i denote the orbital inclination angle with respect to the line of sight. 



A. Mass Quadrupole 



By direct calculations, we obtain the plus mode of quadrupole waves as 



r x h 



Q 



—2 m p a 2 p u) 2 (l + cos 2 i) cos 2(ut + 9 p ), 



(10) 



6 



where the subscript Q denotes a mass quadrupolar part. By using Eqs. (1A4j) - (1A6j) and 
flAT0|) - flAT2l) . Eq. ([I0D is rewritten as 



r x h 



Q 



-m to ta 2 u 2 (l + cos 2 i) 



x 



V2 + ^3) — 2z/ 2 z/ 3 J cos 2ujt + \Tiv\{v2 — ^3) sin 2u;i 



(11) 



In a similar manner, we obtain the cross-mode as 



r x hi 



-4 m p a 2 p uo 2 cos i sin 2 (cut + 9 p ) 



—2mtota 2 ^ 2 cosi 



^1(^2 + ^3) _ 2z/ 2 z/3^ sin 2^ — -\/3^i(^2 — ^3) cos2a;t 



(12) 



For a three-body system in Lagrange's orbit, the frequency sweep due to gravitational 



radiation reaction has been recently obtained as 



io| 



1 d/cw 96 



few dt 
where we define a chirp mass as 



— 7T 8/3 M 5/3 f 8/3 
g 71 J w chirp ./GW' 



(13) 



M, 



chirp 



m t0 



^P P V mtot 



2/3 



^p¥=<1 p q \mtot J 



2/3 / M off x 2/3 



sin 



p P V mt 



a \ 2/3 



3/5 



,(14) 



and Mp denotes the effective one-body mass for which the equation of motion becomes 24] 

d 2 X n 



dt 2 



MfXp 
" I^pI 3 ' 



(15) 



For instance, Mf s is defined as 



off 



(m 2 + m 2 m 3 + m 2 ) 3//2 



(16) 



tot 



By cyclic permutations, M| ff and M| ff are defined. It is worthwhile to mention that this 
frequency evolution equation is the same as that for a binary system [lo| including the 
numerical coefficient 96/5. 
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The chirp mass is expressed in terms of both each mass m p and its initial position angle 
dp. It is rewritten as a function only of the mass ratios by deleting the initial position angles 
as (See Appendix for detailed calculations) 

M chirp = m tot x F(u u u 2 , vz), (17) 

where we define F as 

\v\{v 2 - u 3 ) 2 + z/f(z/ 3 - vi) 2 + vlfa - ^ 2 ) 2 V /S 



(18) 

2 v x v 2 + v 2 v?, + v 3 i>i J 

By using this expression of the three-body chirp mass for Eqs. (TTTjl and (fl2l) . we obtain 



where we define \1/q as 



for 



x \Jv x v 2 + v 2 v z + v z v x cos(2co>t - ^q), (19) 



tan* Q =^, (20) 



Aq = v x {v 2 + v z ) - 2u 2 u 3 , (21) 
B Q = VSvifa - uz). (22) 



We rewrite the cross mode as 



rxh* = -2m^M^y/s C0S i 



x \Iv\Vi + v 2 v% + U3U1 sin(2cjt - \I/q). (23) 

B. Mass Octupole 

We obtain, by direct calculations, the plus mode of quadrupole waves as 

r x ^oct = - Y^Yl m v a l^ sin * 
p 

x [27(1 + cos 2 i) cos 3(ut + 9 P ) + (1 - 3 cos 2 z) cos(c;£ + 6 P )] , (24) 
where the subscript Oct denotes a mass octupolar part. 
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By using Eqs. (R7jl - (R9|l and (lAT3l) - (lAT5|) . one can rewrite Eq. §H§ as 

1 o 



r >< ^Oct 



12 



27(1 + cos 2 i) (3 3//2 z/!Z/ 2 z/3 cos 3ut + (u x — v 2 ) [y 2 — z/ 3 ) (u 3 — u{) sin 3o>t) 



2 .w \/3 



+ (1 - 3 cos i) [ — ^{^2(^2 - Vi) + ^3(^3 - vi)} cosut 



--(^2 - ^){(^i - v 2 ){yi - u 3 ) - 3z/ 2 z/ 3 }sinu;t 



where we use Eq. flH]). 

We obtain the cross mode of mass octupole waves as 



r x hoa 



— m p apU) 3 sin 2i 



x [27 sin 3 (cut + 6 P ) - sin(u;t + 6 P )} . 



(25) 



(26) 



By using Eqs. (TXijl - fjM]) . (|A7l) - ([A9]) and flM3l - flAl~5]) . one can rewrite Eq. ([26]) as 



r x h Oct 



— m?„ t u; sin 2i 
12 tot 



x 



2J^S s / 2 p 1 i/ 2 i/ 3 sin3u;t — (y x — v 2 ) [y 2 ~ ^3) (^3 — ^1) cos3u;t 



y/3 

- ( ^-^{^2(^2 - v x ) + ^3(^3 - ^1)} sin cut 



+ -(z/ 2 - 1/3) {(^1 - u 2 )(ui - V3) ~ 3v 2 v 3 } cosut 



where we use Eq. 



C. Current Quadrupole 



(27) 



Here, we consider current quadrupolar waves. The plus mode becomes 
rxh c = 3 m tot^sm« 



x 



\/3 

— v\{v 2 {y 2 - v x ) + v s (v 3 - v x )} cosut 



--(z/ 2 - i/ 3 ){(i/i - z^X^i - ^3) - 3z/ 2 z/ 3 }sin^t 



(28) 



where the subscript C denotes a current quadrupolar part. 



The cross mode becomes 



rxh* = -m 2 tot uj sin 2i 



x 



\/3 

— v x {v 2 {v 2 - vi) + v 3 (is 3 - z/i)}sincut 



+ -(^2 - - ^2X^1 - ^3) - 3v2V 3 }cosa>* 



(29) 



Both the mass octupolar and current quadrupolar parts are proportional to m t 2 ot u. Hence 
they can be combined as 

1 , . . 



r X ^Oct+C 



and 



-~m tot uj smi 



9(1 + cos 2 i) (3 3//2 z/iz/ 2 z/3 cos 3ut + {y\ — v 2 ) [y 2 — ^3) (^3 — vi) sin 3u;t) 



'\/3 

-(5 + cos 2 i) [ — v x {v 2 {y 2 - vi) + ^3(^3 - ^1)} cos cut 



--(u 2 - ^3){(^i - ^2)^1 - ^3) - 3z/ 2 ^3}sinu;t 



(30) 



r x h 



Oct+C 



9^3 3 ^ 2 z/iz/2^3 sin3u;t — (z/i — v 2 ) [y 2 — ^3) (1^3 — v\) cos3u;t 
— ^{^2(^2 - ^1) + ^3(^3 - ^1)} sin cut 



+ 2(^2 - ^){(^i - v-i){y\ - ^3) - 3z/2^ 3 }coscjt 



(31) 



It is natural that one can recover the wave forms to a binary system from Eqs. (|30|) and 
(l3Tj) when a third mass vanishes, say 1/3 = 0. 

Figure [2] shows wave forms due to mass quadrupole, octupole and current quadrupole 
parts that are expressed above. 



IV. PARAMETER DETERMINATIONS 



The frequency 00 can be directly determined as uo = 7t/gw from measurements of the 
frequency of the mass quadrupolar wave few- What we have to do is to determine the 
remaining six quantities as m tot , v>i, u 2 , z/ 3 , r and i. We have an identity v\ + v 2 + u 3 = 1. 
Hence, we need to find out five more relations for parameter determinations. 
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FIG. 2: Gravitational waves for a three-body system in Lagrange's orbit. Log-dashed (blue), 
short-dashed (red) and dotted (green) curves denote mass quadrupolar, octupolar and current 
quadrupolar parts, respectively. The sold (black) one denotes the total wave forms. The vertical 
axis is in arbitrary units and time in the horizontal axis is normalized by the orbital period. For 
simplicity, the mass ratio is assumed as mi : rri2 : m% = 1:2:3, though stability arguments prefer 
much larger ratios 



24l ]. In order to exaggerate differences between each component of waves, we 
assume a mildly relativistic case of a = 100m to t) which corresponds to v/c ~ 0.1. The mass 
octupolar part makes a relatively large contribution mostly because of including a large numerical 
coefficient (as 27 in Eq. (|25p ). whereas the current quadrupolar part is much smaller and changes 
slowly with time because it has no 3a; but only the u part. 

The frequency of mass quadrupolar waves is 2u, whereas that for the combination of 
mass octupolar and current quadrupolar ones is either u or 3u. By using this difference 
of frequency dependences, therefore, one can pick up the quadrupolar waves from observed 
signals. By comparing amplitudes of the + and x modes, we obtain 

Am P Oq) = 2 cos i 
Amp (/iq) 1 + cos 2 % 

where Amp denotes the amplitude of waves. This relation is the same as the well-known 
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one for binaries. The L.H.S. of Eq. (1321) can be determined by observations and thus the 
R.H.S. tells us the inclination angle i. 

Through observing the frequency sweep of the mass quadrupolar part (See Eq. (JT31 ). 
one can determine the chirp mass as 

M ^=[^]m J jr) • (33) 

where /gw = 2/ for / = u/2ti. 

Next, let the amplitude of each wave component observed separately. From Eq. (|19[) . 
first, we obtain that for the mass quadrupolar part as 

2 



r x 



X^l^ + V 2 V 3 + ^3^1, (34) 



which can be solved for the distance r as 



Amp (hX) 

For parameter determinations, however, it is useful to rewrite this as 
^totV^i^ + ^2^3 + ^3^1 _ Am P ( h o) 



r 2M^// 3 (l + cos 2 : 

where the R.H.S. is known up to this point. 

By using Eq. (|30l) . the magnitude of ft-Qct+c a ^ 3o> becomes 

9 

r x Amp (/*oct+cU) = ^m^wsin^l + cos 2 i) 



(35) 



(36) 



x^27z/ 2 z/ 2 2 z/f + {v x - v 2 )\v 2 - v 3 )\v 3 - v x )\ (37) 

where the subscript of 3a; in the L.H.S. means a 3a; part. This equation is rewritten as 

m\ ^21v\vlvi + {v x - v 2 f{v 2 - v 3 f{v 3 - v x f 



r 



Amp (h+ ct+c \ 3u} ) , (38) 



9 oj sini(l + cos 2 i 
where the R.H.S. is known up to this point. 

By using Eq. (|30l) . the magnitude of fiQ Ct+c at to becomes 

r x Amp Ooct+cU) = \ m lo^ sin *( 5 + c °s 2 

x [3{i/i^ 2 (vi - v 2 ) + v 2 v 3 {v 2 - u 3 ) + v 3 v x {v 3 - i>i)} 2 
+ {vx-V2)\v2-vi)\v?,-vx) 2 \ Xl \ (39) 
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where the subscript of u in the L.H.S. means a u part. This equation is rewritten as 



m 



tot 



X 



3( v 1 v 2 [u 1 - v 2 ) + v 2 v 3 {y 2 - vg) + v&iiyz - Vx. 

1/2 



+ (^1 - V 2 ) {V 2 - I/3) (^3 - ^l)' 
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where the R.H.S. is known up to this point. 
Equation (130]) is rewritten as 



(40) 



r X ^Oct+C 



9(1 + cos 2 i)\jHvlvlvl + ( Vl - v 2 f[v 2 - u 3 ) 2 {u 3 - v x f cos(3wt - ^) 
- -(5 + cos 2 i) ^3{z/iz/ 2 (^i - v 2 ) + ^3(^2 - u 3 ) + vzv x {v 3 - z/i)} 2 



+ {V X - V 2 ) 2 (V 2 - I^) 2 (z/ 3 - ^l) 2 



1/2 



cos (cut — ty u 



tan 



where phases of /iQct+c a ^ w anc ^ 3w are defined as 

(u 2 - vz){{vi - v 2 ){vi - z/ 3 ) ~ 3z/ 2 ^ 3 } 
y/3vi{v 2 (v 2 - ui) + 1/3(1/3 - ui)} 

{v\ - ^2) (^2 - ^3) (^3 - vi) 



and 



tan* 



3u> 



(41) 



(42) 



(43) 



3 3 / 2 z/iZ/ 2 z/ 3 

respectively. These phases are not directly observable. Time lags between the mass 
quadrupole wave with frequency 2uj and the combination of mass octupole and current 
quadrupole waves with to (or 3u) can be measured. They are defined as 



At, 



At 



3u> 



t u — tq, 
2u 

2tf„ - 3*r 



6co 



respectively, where £q, t u and t 3w are defined as 



+ — Q 

tQ = 2Z7' 
t = — 

J 3w = t; — ■ 

6UJ 



(44) 

(45) 

(46) 
(47) 
(48) 
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It is worthwhile to mention that At^ and At^ are observable and thus gauge- invariant, 
whereas £q, t w and t ?jU) are gauge-dependent in a sense that they rely on a degree of freedom 
for choosing an initial time "to" • 

In principle, time lags in the cross mode, phase of which is different from that of the 
plus mode by 45°, are the same as those in the plus mode, and thus bring no additional 
information on parameter determinations, though they may play a supplementary role in 
improving accuracy of a practical data analysis. 

As a result, one can determine five quantities m tot , u%, u 2 , v% and r from Eqs. ([7j), (jTTl) . 
( 1361) . ( 1381) . ( f40|) . ( |4~4"1) and ( 145]) . In principle, five out of the seven equations are sufficient 
for parameter determinations if the source is known as the Lagrange's solution a priori. It 
should be noted that two remaining equations are never redundant but play an important 
role in checking whether a source is the particular three-body system or not. If the remaining 
equations are satisfied by the determined parameter values, one can safely say that the source 
is in Lagrange's orbit. If not, it could be other systems. 

Figure [3] shows a flow chart of the parameter determinations and possible source tests 
that are discussed above. 

V. CONCLUSION 

We have considered the three-body wave forms at the mass quadrupole, octupole and 
current quadrupole orders, especially in an analytical method. By using the derived ex- 
pressions, we have solved a gravitational-wave inverse problem of determining the source 
parameters to the particular configuration (three masses, a distance of the source to an 
observer, and the orbital inclination angle to the line of sight) through observations of the 
gravitational wave forms alone. We have discussed also whether and how a binary source can 
be distinguished from a three-body system in Lagrange's orbit or others, and thus proposed 
a binary source test. 

To be more precise, we should take account of post-Newtonian corrections to both wave 
generation and propagation beyond the linearized theory. It is interesting to consider dif- 
ferent shapes of orbits and different numbers of bodies (N = 4 or more). This is a topic of 
future study. 
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TABLE I: List of quantities characterizing a system in this paper. 



Symbol Definition 

T Orbital period 

uj Angular velocity (=2ir/T) 

a Edge length of a Lagrange's equilateral triangle 

a p Distance of each body from their center of mass 

6 p Initial angular position of each body 

m p Mass of each body 

up Mass ratio of each body 

TOtot Total mass 

-^chirp Chirp mass 

F Ratio as M^pin^ 

r Distance of a source from an observer 

/iq' x + or x mode of mass quadrupolar waves 

Hq'^ + or x mode of mass octupolar waves 

Hq' x + or x mode of current quadrupolar waves 

h + ' x h + ' x +h + ' x 

n Oct+C n Oct + n C 

Amp Amplitude of a wave component 

\£q Phase of mass quadrupolar waves with 2lo 

^ Phase of waves with uj 

^Zu) Phase of waves with 3uj 

tQ Time corresponding to the phase 

t u Time corresponding to the phase ^ u 

Time corresponding to the phase 

At^ Observable time lag defined as t w — tQ 

Atzu Observable time lag defined as t% w — tQ 
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APPENDIX A: LAGRANGE'S SOLUTION IN COMPLEX VARIABLES 
1. Correspondence between position angles and mass ratios 

Let each mass located at each vertex of a regular triangle with the side length denoted 
as a. For simplicity, the center of the complex plane z = x + iy is chosen tentatively 
as that of the triangle. Without loss of generality, we assume initial positions of three 
masses as a, ae i27r//3 , ae t4n ^ 3 . The center of mass for the three masses is located at zcm = 
a {y x + u 2 e i2w / 3 + z/ 3 e i47r / 3 ) . 

When we consider gravitational waves, it is convenient to choose the coordinates origin 
as the mass center so that the dipole moment can vanish. In order to do so, we make a 
translation as z — > z — zcm, for which the initial position of each mass becomes 

zn = ^-a (y3(u 2 + u 3 ) ~ (v 2 - v 3 )ij , (Al) 
= -^j-a (VSvi - {yx + 2iA,)i) , (A2) 



2 



> s = ~<i (V3vi + (i/i + 2v 2 )i) , (A3) 



where the subscript I means the initial values. The magnitude of zn, z\ 2 and z I3 , which are 
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nothing but three masses' orbital radii around their center of mass, become 



a i = \Zll\ 



= ay / 3(z/ 2 2 + z/| + z/ 2 z/ 3 ), (A4) 



a 2 = \z\2\ 



= ay / 3(z/ 3 2 + z/ 1 2 + z/ 3 z/i), (A5) 

«3 = |^I3 1 

= ay / 3(z/ 1 2 + ~v\ + v x v 2 ), (A6) 

which respect symmetry for cyclic changes asl— >2— >3— >1. 

We obtain a relation between an position angle 6\ (with respect to the center of mass) 
and the mass ratios as 



e L = 



\zn\ 

V3{v 2 + ^3) - ("2 - vz)i 

1\Jv\ + Z/f + V 2 Vi 



(A7) 



Similarly, we obtain 



e*» = _ ^~^ + ^)\ (A8) 
2a/z/| + z/ 2 + 1/3^1 



2a/^ 2 + i/f + z^ 2 

By using the complex representations, double-angle and triple-angle relations, which are 
useful for later computations, are obtained as 



2(i/f + i/f + v 2 v z ) 

2,0., ("I ~ - 2z/ 3 z/i) - y^3i/i(i /i + 2u 3 )i 
2(z/| + z/ 2 + v z v x ) 

2i e 3 = K - 2v l - 2 "i^) + v^i^i + 2v2)i (Am 

2(z/ 2 + z/ 2 + v x v 2 ) ' 1 ' 



e2 i8! = \ u 2 Z "3 I ^^2^3; - V o\T2 ~ ^3^ (A10) 



e 2i(? 2 = ^1 ~ ^3 ~ ^3^1 j ~ V T ^3)0 (All) 



and 



,3i0i 



3 3/2 z/ 2 1/3(1/2 + 1/3) - (z^2 - ^3) (^2 + 2z/ 3 )(2z/ 2 + z/ 3 )i 



2( 


A + »i 


+ V 2 Vzfl 2 


3 3 / 2 z/ 3 z/i(z/ 3 + z/i; 


! + (i/l - 


- vz)iy\ + 2i/ 3 )(2i/i + z/ 3 )i 


2( 


;^ 3 2 + ^ 2 


+ ^l) 3/2 


3 3 / 2 z/!Z/ 2 (z/i + v 2i 


1 - (y\ - 


- ^2X^1 + 2i/ 2 )(2z/ 1 + i/ 2 )i 



(A13) 



C o/..2 , ..2 , .. .. \3/2 ' V^ 1 ^) 

^ = ^(z/f + ^+^W) 3 / 2 ~' (A15) 
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respectively. 

Some useful relations are obtained as 



A c = } ] v P \zi p \ s cos d p 
p 

9 



-vM + vl -vi(v2 + v 3 )]a 6 , (A16) 

= v P \zip\ 3 sin Op 
v 

33/2 



V2 ~ v$[v x {y x - v 2 - vz) - 2v 2 v 3 ]a 3 , (A17) 



2 

where we use Eqs. (IA7p -f lA~9]) . Therefore, we obtain 



^ v P \zip\ 3 cos(ut + 6 p ) = JA C + Alcos(ut + a), (A18) 
p 

where we define 

a = arctan ( -7— ] • (A19) 



.4 



Next we consider 3u parts. We obtain 

B c = y^z/ p |2:ip| 3 cos 39 p 
p 

= 27z/!Z/ 2 z/ 3 a 3 , (A20) 

B s = v p \zi p \ 3 sin 39 p 
p 

= -3 3/2 (z/! - v 2 ){v 2 - z/ 3 )(z/ 3 - v x )a\ (A21) 
where we use Eqs. (IA13[) -f lAT5|) . Therefore, we obtain 



^^|zi p | 3 cos3(c;t + ^) = JB C + 5§cos(3u;t + /3), (A22) 
p 

where we define 

(3 = arctan (^-) . (A23) 



c 



2. Chirp mass 

By straightforward but lengthy calculations, one can show several identities as 

2/3 



M eff \ 

5^ v p ( = v ^ + ^2^3 + Wi, (A24) 

^ 1 m tot ' 



) 
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2/3 

) 



■ p U q -Y]Vp[ = V\V2 + + V3V1, (A25) 



p^q p 

2 / 3 / wrfF \ 2/3 



2 L¥g — — sin2 (^ = 3 ^i^3, (A26) 

p _, g V mto v V mtot / 

where Eq. is frequently used. 

By substituting these relations into Eq. ffT4l) . we obtain 

= i[z/ 2 (z/ 2 - z/ 3 ) 2 + z/ 2 > 3 - ^i) 2 + ^ 3 >i - ^) 2 ]- (A27) 

Therefore, we can prove that the chirp mass defined by Eq. (fl4l) is expressed in terms of 
only the mass ratios as Eq. (FT?]) . 
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FIG. 3: Flow chart of the parameter determinations and source tests. By using equations derived 
in this paper, the source parameters can be determined through gravitational-wave observations 
alone. In addition, a binary source can be distinguished from a three-body system in Lagrange's 
orbit or others. 
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